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Abstract. The optical-ultraviolet continuum from active 
galactic nuclei (AGN) seems to originate from optically 
thick and/or thin disks, and occasionally from associated 
circumnuclear starburst regions. These different possible 
origins can, in principle, be discriminated by observations 
of gravitational microlensing events. We performed nu- 
merical simulations of microlensing of an AGN disk by a 
single lensing star passing in front of the AGN. Calculated 
spectral variations and light curves show distinct behavior, 
depending on the nature of the emitting region; time vari- 
ation of a few months with strong wavelength dependence 
is expected in the case of an optically thick disk (standard 
disk), while an optically thin disk (advection- dominated 
disk) will produce shorter, nearly wavelength- independent 
variation. In the case of an associated circumnuclear star- 
burst region much slower variations (over a year) will be 
superposed on the shorter variations caused by microlens- 
ing of the disk. 

Key words: accretion, accretion disks — active galactic 
nuclei — microlensing 



1. Introduction 

There has been long discussion regarding the central struc- 
ture of active galactic nuclei (AGN). Generally, it is be- 
lieved that there is a supermassive black hole and a sur- 
rounding accretion disk. However, it is quite difficult to 
obtain direct information about the center of AGN, the 
accretion disk, because the disk size is far too small to re- 
solve. As for cosmological objects, e.g., quasars, the (an- 
gular diameter) distance is of the order of ~ IGpc ~ 
3 x 10 27 cm. If a black hole placed at the center has a 
mass of ~ 10 8 Af Q , the disk size (say ~ 1000r g , where r g 
is Schwarzschild radius) will be ~ 3 x I0 16 cm. Hence, 
apparent disk size (#disk) is only ~ I0 _11 radian ~ l//as. 
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What is usually done is to observe the spectra, which 
are the sum of local spectra over the entire disk, and to 
make a fitting to the observed spectra by theoretically in- 
tegrated spectra. However, there remain ambiguities as to 
the spatial emissivity distribution. The standard picture is 
that the so-called UV bump is attributed to blackbody ra- 
diation from the standard- type, optically-thick disks (i.e., 
Shakura & Sunyaev 1973). Such a fundamental belief may 
be reconsidered by the recent HST observations that show 
no big blue bumps (e.g., Zheng et al. 1997). 

Instead, microlensing can be used as a 'gravitational 
telescope' to resolve the disk structure (see Blandford & 
Hogg 1995; originally proposed by Chang & Refsdal 1984), 
since a typical Einstein-ring radius of a stellar mass lens is 
very small; e.g., 6*e ~ I0~ n radian for cosmological distant 
sources. With such a gravitational telescope we may be 
able to obtain information about the disk emergent spec- 
tra as a function of the distance from the central black 
hole. 

When a point source object is lensed, its luminosity 
is amplified without any color change. This feature is 
called 'achromaticity' and is discussed in many literature, 
initiated by Paczyhski (1986; for a review see Narayan 
& Bartelmann 1996). But, in some actual situations, we 
cannot treat a source as a point and the effect of a fi- 
nite source size should be taken into account. Grieger, 
Kayser & Refsdal (1988) examined such finite size ef- 
fects. They simulated microlensing light curves for two 
different source profiles: stellar like and accretion disk 
like. The light curves show distinct shapes from those of 
point source calculations. However, they only considered 
monochromatic intensity variation. Later, Wambsganss & 
Paczyhski (1991) pointed out that if the finite size source 
has inhomogeneous structure so that each part produces 
a different spectrum, lensing light curves have wavelength 
dependence; i.e., chromatic features in microlensing events 
should appear. This chromatic effect has been investi- 
gated by several authors who calculated microlensing light 
curves of AGN accretion disk based on simple disk models 



2 



A. Yonehara et al.: Microlens diagnostics of accretion disks in active galactic nuclei 



(e.g., Rauch & Blandford 1991; Jaroszyriski, Wambsganss 
& Paczyriski 1992). 

The purpose of this paper is to obtain some criteria 
for distinguishing disk physics based on more realistic disk 
models, we adopt the optically-thick, standard-type disk 
and the most successful, optically thin disk model, the 
advection-dominated accretion flow (ADAF; Narayan & 
Yi 1995; Abramowicz et al. 1995). So far nobody (except 
Yonehara et al. 1998) has yet calculated X-ray or radio 
emission in the microlens events, due probably to the lack 
of reliable disk models. However, such calculations arc of 
great importance since X-ray and radio emission seems 
to originate from the vicinity of a putative black hole. 
In section 2 we describe our methods of simulation. In 
section 3 our results are displayed. Section 4 is devoted to 
discussion. 

2. Method of the simulation 

2.1. Accretion disk models 

First, we describe the adopted accretion disk models. We 
consider two representative types of disks: 
1. Optically thick disk: 

The first case we consider is the geometrically-thin and 
optically-thick standard disk (Shakura & Sunyaev 1973). 
Since we are concerned with disk structure on scales rang- 
ing from a few tens to a few thousands of r g , we ignore 
relativistic effects, such as the gravitational redshift by the 
black hole, for simplicity. The temperature distribution of 
the optically thick standard disk is then given by (Shakura 
& Sunyaev 1973) 



T(r) = 2.2 x 10 5 
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where M is the mass accretion rate and M is the black 
hole mass. The inner edge of the disk, n n , is set to be 
ri n = 3r g , the radius of the marginally stable last circular 
orbit around a non-rotating black hole. 

Since the disk is optically thick, we can assume black- 
body radiation. Inserting the temperature profile T(r) into 



AL{v,r) = 2 • AttB v [T(r)] AS* cos a, 



(2) 



we can calculate the luminosity of a part of the disk with 
the surface element, AS = Ar ■ rA(p, at a radius r and 
a frequency where B v is a Planck function (e.g., see 
eq.(1.51) of Rybicki & Lightman 1979), and a is the incli- 
nation angle of the disk. For simplicity, here, we assume 
the disk is face-on (a = 0). It might be noted that the 
radiation spectra do not depend on the viscosity in this 
standard-type disks, since the total emissivity and the ef- 
fective temperature are solely determined by the energy 



balance between the radiative cooling and viscous heating 
(which originates from release of gravitational energy). 
2. Optically thin disk 

The second model is the optically thin version of the ad- 
vection dominated accretion flow (ADAF). We used the 
spectrum calculated by Manmoto, Mincshigc, & Kusunose 
(1997) and derive luminosity AL(v 1 r) in the same way as 
we did for an optically thick disk: 



AL(u, r) = 2 • e v (r)AS, 



(3) 



where e v is surface emissivity which include the effect of in- 
clination angle (details are given in Manmoto et al. 1997). 
In this paper, we also set a = (face-on) same as the case 
of optically thick disk. Here, included are synchrotron, 
bremsstrahlung, and inverse Compton scattering of soft 
photons created by the former two processes (see Narayan 
& Yi 1995 for more details). Optical flux is mainly due to 
Comptonization of synchrotron photons. 

In this sort of disks (or flows), radiative cooling is in- 
efficient because of very low density. As a result, accreting 
matter falls into a central object, hardly losing its internal 
energy (which is converted from its gravitational energy 
through the action of viscosity) in a form of radiation. 
Hence, although the total disk luminosity is less than that 
of the standard one for the same mass-flow rate, the disk 
can be significantly hotter, with electron temperatures be- 
ing of the order of 10 9 K or more, thus producing high 
energy (X-7 ray) photons. Photon energy is thus widely 
spread over large frequency ranges. Since emission inside 
the last circular orbit (at 3r g ) is not totally negligible in 
this case, we solve the flow structure until the event hori- 
zon (r g ) passing through the transonic point, and con- 
sider radiation from the entire region outside the horizon, 
although the contribution from the innermost part is not 
dominant. Unlike the optically thick case, the viscosity ex- 
plicitly affects the emissivity in this case, since radiative 
cooling depends on density and is no longer balanced with 
viscous heating (and with gravitational energy release). In 
the present study, we assign the viscosity parameter to be 
a = 0.1. 

2.2. Microlensing events 

Although single lens approximation is not always appro- 
priate to the microlensing events, we, here, use this simple 
method because of making the situation clear (for extend 
source effect of microlensing by single lens object, e.g., 
Bontz 1979). 

Generally, the microlensed images are split into two (or 
more) images. When the separation between these images 
is too small to resolve, however, the total magnitudes of 
all the split images vary due to this microlensing caused 
by the relative motion of lens and source. 

The properties of microlensing have been extensively 
examined for simple cases where both the lens and source 
are regarded as being points (e.g., Paczyhski 1986). In our 
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simulation, we include an effect of occultation by a lens ob- 
ject which is recently included the simulation of microlense 
light curves by Bromley (1996) for reality. In other words, 
we, here, consider that the radius of the lens object (i?i en s) 
is same as the sun, i.e., i?i cn s = ^0 = 6.96 x 10 10 cm, and 
calculate image positions of the all cell. If the image posi- 
tion of a cell from the lens object is smaller than the radius 
of lens object, we exclude a contribution from that image. 
Therefore, the total amplification (or magnification) fac- 
tor A for microlensing is exactly expressed as 



Lens Object 



lens 
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(otherwise), 



where u corresponds to the angular separation between 
lens and source in the unit of the Einstein-ring radius 
and always greater than or equal to zero, i.e., u > 0, and 
#imagc,± shows that the image positions of two microlensed 
images from a lens object in the unit of radian, i.e., 



1 
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This function as shown in cq. ^ is a monotonically increas- 
ing function of u; A(u) ~ 1.34, ~ 10, ~ 100 for u = 1.0, 
0.1, and 0.01. 

Here, the Einstein-ring radius is 
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where Mi cns is the mass of a lensing star, c is the speed 
of light, G is the gravitational constant, and D\ s , D os , and 
D a \ denote the angular diameter distances from lens to 
source, from observer to source, and from observer to lens, 
respectively (Paczyhski 1986). Adopting an appropriate 
cosmological model, we can calculate these angular diame- 
ter distances in terms of redshifts z\ s (from lens to source) , 
z os (from observer to source), and z \ (from observer to 
lens). In the present study, we assume the Einstein-de Sit- 
ter Universe, in which we have 



2c 1 
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(l + z x ) 1/2 
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(e.g., Padmanabhan 1993), where Hq is Hubble's constant, 
and the subscript 'x' stands for 'Is', 'os', or 'ol\ 

Since we now consider microlensing effects on an ex- 
tended source, we must integrate the magnification effects 
over the entire source plane. We thus, as a next step, need 
to calculate the angular separation u between a part of 
the source in question and the lens center (see figure 1), 
from which the amplification factor, A(u), can be found. 
In the present study, we first divide an accretion disk plane 
into azimuthal (</>) and radial (r) elements. The azimuthal 
coordinate is equally divided into ~ 1000 segments (</>j), 




Fig. 1. Schematic view of calculation for the separation 
between each cell and lens. 

while the logarithmically scaled radial coordinate is also 
equally divided into ~ 1000 segments (ri). We then calcu- 
late the observed flux at frequency v from a cell between 
and ri + \ in the radial direction and between cf>j and 
in the azimuthal direction. Noting that photons emitted 
by the source at a frequency v are observed at the fre- 
quency z//(l + z os ), we write the observed flux as 



AF obs (^;rj,^) ~ A(u)AF ohStQ [v(l + z os ),r l; ^ 



(8) 



where A{u) is the amplification factor given by equation 
(4). Moreover, u is explicitly given by 



1 cos <b\ ens 
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where a is the inclination angle of an accretion disk (in 
this paper, as indicated in §2.1, we set a — 0), and 6^ and 
</>icns represent the radial and the azimuthal positions of 
the lens in units of radians, respectively. Schematic view 
of the calculated disk is shown in figure [l]. The intrinsic 
flux in the absence of a microlensing event, Ai 7, bs,o( !/ j r )i 
is calculated from AL(v, r) according to 



Afobs,o(f, r) = 



AL(v,r) 



47rL>2 s (l + z 



(10) 



and AL(v,r) depends on the disk model (see eq. (Q) and 
(||)). In the present study, since we neglected the rela- 
tivistic effect (e.g., beaming, gravitational redshift etc.), 
AL(i>; r, <fi) has no azimuthal dependence. 

Summing up AF ^ s {v; r, cj)) over the entire disk plane 
from the inner boundary (n n ) to the outer boundary 
('"out)! we obtain the total observed flux at frequency u, 
or the spectrum of the microlensed accretion disk. 
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3. Results of the calculations 

3.1. Model parameters 

For clarity, we consider a lensed quasar, Q 2237+0305 (the 
Einstein cross, e.g., Huchra et al. 1985). Irwin et al. (1989) 
reported an increase of the apparent luminosity of image 
A by ~ 0.5mag on a timescale of a few months, which was 
nicely reproduced by a model of microlensing (Wambs- 
ganss, Paczyhski, & Schneider 1990). At present, at least 
five microlens events have been reported so far on the split 
images (Irwin et al. 1989; Corrigan et al. 1991; Houde 
& Racine 1994). Although the single lens treatment for 
Q 2237+0305 may be a poor approximation (see discus- 
sion), we just take this object as one test case and derive 
quite general conclusions regarding the generic observa- 
tional features of microlensed accretion disks. 

This object is macrolensed by a foreground galaxy, 
and the source is split into four (or five) images. Now, 
we suppose that one of the macrolensed image is further 
microlensed by a star in the foreground galaxy causing 
the macrolensing. Although the macrolensed images are 
actually affected by two effects, so called, 'convergence' 
and 'shear', we here neglect both for simplicity and thus 
assume that the macrolensed image is neither amplified 
nor distorted. The redshift parameters are 



1.675, z a \ = 0.039 



(11) 



(cf. Irwin et al. 1989), and substitute them into the rela- 
tion 



1 



1 



1 + Zo\ 



(12) 



which can be obtained from flux conservation of radiation 
(e.g., Padmanabhan 1993), we get 



zis = 1-575. 



(13) 



As for Hubble's constant, we adopt the value obtained by 
Kundic et al.(1997), H Q ~ 60km s^Mpc" 1 . The adopted 
disk parameters are the outer edge of the disk, r out = 
10 3 r g , the mass of the central black hole, M — lO 8 Af0, 



10 26 g s 1 for the stan- 



and the mass-flow rates, M = 
dard disk and M = 10 22 g s _1 for the optically thin disk. 
We take different mass-flow rates for two cases because 
of calculation convenience. This is not, fortunately, a se- 
rious problem, since it has been demonstrated that the 
spectral shape of the ADAF is not very sensitive to M, 
although the overall flux level does change. We will not 
discuss the absolute flux in the following discussion, and 
focus on the relative spectral shape and relative flux vari- 
ation. For other basic parameters characterizing the disks, 
see §2.1. 

For simplicity, we assume that the accretion disk is 
face-on; i.e. i = 0, to the observer. Actually, based upon a 
grand unified paradigm of AGNs, quasars are not far from 
face-on views. We also assume the mass of a lensing object 
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Fig. 2. Microlensed spectrum of a standard accretion disk 
(upper) and an optically thin accretion flow (lower). In 
each panel, the microlensed spectra are plotted for angular 
separations between the disk center and the lensing star 
to be u = 1.0 (dotted line), 0.1 (dashed line), and 0.01 
(long dashed line), together with the one in the absence 
of microlensing (solid line). The lensing mass is Mi ens = 
1.0M Q . 



to be Miens = I.OMq or O.IMq. The remaining variable 
is the angular separation (it) between the lens and a part 
of the source in question. 

Adopting these parameters, apparent angular accre- 
tion disk size (#disk) and Einstein ring radius (#e) is 
6> disk - 6.8 x 10~ 12 rad and E - 3.3 x 10 _n rad. These 
values are comparable and we, here, can use microlensing 
event as effective 'gravitational telescope' (Blandford & 
Hogg 1995). 



3.2. Spectral variation 

First, we calculated the microlensed spectra of each of two 
types of accretion disks for the angular separation between 
the lens and the source of u = 1.0, 0.1, and 0.01, respec- 
tively. The results are shown in figure ^| for the standard 
disk (in the upper panel) and for the optically thin disk 
(in the lower panel), respectively. For the former disk not 
only the disk brightening, but also substantial spectral de- 
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Fig. 3. Contribution to the total integrated spectra from 
different parts of disks for the cases of a standard (upper) 
and optically thin (lower) disk, respectively. Each figure 
displays the contributions from the inner region of r < 
10r g (long dashed line), the intermediate region of 10r g ~ 
100r g (dashed line), and the outer region of 10 2 r g ~ 10 3 r g 
(dotted line). 



Fig. 4. The radial distribution of the emergent flux with 
several wavelengths for the case of a standard-type disk 
(upper) and an optically thin disk (lower), respectively. 
The adopted bands are U-band (3650 A, by the solid line), 
V-band (5500 A, by the dotted line), I-band (8800 A, by 
dashed line), K-band (22000 A, by long dashed line), and 
X-ray (5 keV, by the dot-dashed line). 



formation is produced by microlensing, when the angular 
separation is small (Rauch & Blandford 1991; Jaroszyhski, 
Wambsganss & Paczyhski 1992). The smaller the angular 
separation is, the larger becomes the modification of the 
observed spectrum. In other words, the amplification fac- 
tor depends on the frequency (wavelength) of the emit- 
ted photons. This gives rise to frequency-dependent, mi- 
crolensing light curves (see §3.3). As for the optically thin 
disk, in contrast, there will not be large spectral modifi- 
cation by microlensing, but the flux is amplified over wide 
frequency ranges. 

To understand such spectral properties, we divide the 
disk plane into three parts; r < 10r g , 10r g < r < 100r g , 
and 10 2 r g < r < 10 3 r g , and plot the spectrum from each 
part of the disk, together with the integrated spectra, in 
figure ||. 

In the case of the optically thick accretion disk, pho- 
tons emitted from different radii have different dominant 
frequencies because the local temperature of the disk has 
a strong radial dependence, T oc r~ 3 / 4 (Eq. [1]). Note 



that without lensing the total spectrum has a peak at 
v ~ 10 153 Hz, but its low-frequency part has a smaller 
slope [yFv oc ^ 4 / 3 ) than that of blackbody radiation 
(vF v oc z/ 3 ), thus a shoulder like feature extending to 
v ~ 10 14 Hz (~ 300001) being formed. This is because 
the outer cooler portions also contribute to the spectrum. 
For u = 0.01, therefore, the emission only from the inner, 
hottest part is strongly magnified, much strengthened a 
peak at v ~ 10 15 3 Hz (~ 1500A) and a shoulder like fea- 
ture at v ~ 10 14 Hz (~ 30000A) appears to be weakened. 
Consequently, microlensed spectrum became a form with 
one sharp peak. 

To help understanding such situation, we also plot the 
radial distribution of emitted photons with several wave- 
lengths in figure f§ The upper panel of figure ¥& shows 
rather gradual change of each flux over a wide spatial 
range. This reflects the fact that viscous heating and ra- 
diative cooling are balanced in the standard disk. Since 
the potential energy only gradually changes with the ra- 
dius, so is the viscous heating rate and is the radiative 
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cooling rate. Moreover, the lower a photon frequency is 
(or the longer the wavelength is), the wider becomes the 
parts of the disk which generate photons. 

In the optically-thin accretion disk, on the other hand, 
the optical-UV flux is totally dominated by that from the 
inner region (r < 10r g ). The lower panel of figure^] clearly 
shows that almost 100% of optical flux originates from 
the region inside 10r g , and that the contribution from the 
outer parts is smaller by several orders. This is possible in 
ADAF, since radiative cooling is no longer directly related 
to the shape of the gravitational potential well (which has 
a rather smooth profile). The reason why optical flux is 
produced in a rather restricted region inside 10r g is due to 
enhanced synchrotron emission in the vicinity of the black 
hole, where density and pressure (and thus magnetic pres- 
sure) are likely to be at maximum within the flow. Nearly 
irrespective of the u value, as a consequence, it is always 
the innermost region (within 10 r g ) that dominates the 
entire disk spectrum. The amplification of the total emis- 
sion just reflects that of the emission from the innermost 
parts. Therefore, the overall flux will be amplified by lens- 
ing without large spectral changes. 

Figure |J also shows that X-rays are produced over 
a wider range; the contribution from the outer part at 
100 — 1000r g is not negligible at frequency of / ~ 10 18 (a 
few keV). This is because of a bremsstrahlung emission 
from high-temperature electrons (with T e ~ 10 9 K being 
maintained at r ^ 100r g , see Manmoto et al. 1997). We 
can thus predict that small frequency dependence will be 
seen in X-ray ranges in this model. 

Such distinct spectral behaviors give rise to different 
light curves of the two disk models (see below) . 

3.3. Light curves 

Using the results of the previous subsection, we now cal- 
culate the light curve during the microlensing events (see 
fig. |l|) . We continuously change the angular separation be- 
tween the lens and the center of the accretion disk (9d) 
according to 



Case 1 
(Fig. 6) 



9 d = 6 d (t) = • 



D. 



ol^E 



1/2 



(14) 



and calculate the flux at each frequency, where b is the 
impact parameter (in the unit of 0e), v t is the transverse 
velocity of the lens and t is the time. Note that vt also 
includes the transverse velocity due to the peculiar motion 
of the foreground galaxy relative to the source and the 
observer. We calculated three cases (see fig. ||). 

Figure |^ shows the microlensing light curves of opti- 
cally thick (upper panel) and thin (lower panel) disks for 
the case with (6, M lcns ) = (0.01, 1.0Af Q ) (Case 1). The 
light curve of the optically thin disk is still achromatic 
because of its flat temperature distribution, whereas the 



Case 2 
(Fig. 7) 




Case 3 
(Fig. 8) 



1000 r g /D os 6e 




Fig. 5. Schematic view of the microlensing events which 
we have calculated. These figures show the relative sizes 
of the accretion disk, impact parameter, tracks which lens 
object passes, and the Einstein-ring radius. In each panel, 
shaded regions correspond to the disk with a radius of 
100r g . Top panel (Case 1) shows microlensing event of 
(b, Mi ens ) = (0.01, I.OMq). Left, thick outer and inner cir- 
cles corresponds to the radius of 100r g and 10r g , respec- 
tively. A right circle corresponds to ~ 0.01#e- Inside this 
circle, the amplification factor caused by microlensing is 
more than 100. Middle panel (Case 2) shows microlensing 
event of (b, M lcns ) = [0.01(M Q /Mi en8 ) 1/2 , O.1M )]. The 
types of left circles are the same as those in the top 
panel. Bottom panel (Case 3) shows microlensing event 
of (6, Miens) = (0.1,1.0M Q ). Left, thick outer and inner 
circles corresponds to 1000r g and 100r g of accretion disk. 
A right circle corresponds to ~ 0.1#e, amplification factor 
caused by microlensing is more than ~ 10. 



light curve of the optically thick one shows strong fre- 
quency dependence. 

Figure^ is again useful to understand this chromaticity 
and achromaticity. Since flux distribution is similar among 
different optical fluxes in the case of optically thin disks, 
each optical flux is amplified in a similar way and thus 
the color does not change; i.e. achromatic feature appears. 
X-ray emitting region is wider than that of optical flux, 
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Fig. 6. The microlensed light curves of the standard disk 
(upper) and for the optically thin flow (lower) . The param- 
eters are (b, M lcns ) = (0.01, 1.0M Q ) (Case 1). The different 
type of lines correspond to the same as in figure ||. 

small chromaticity may appear when we compare optical 
and X-ray light variations. 

In an optically thick disk, conversely, low-energy pho- 
tons are created in a wider spatial range, compared with 
high-energy photons. If the microlensing event occurs, 
therefore, low-energy emission will be the first to be am- 
plified, followed by the significant amplification of high- 
energy emission. Note that each optical flux has the same 
level at the beginning (at —120 day) in the upper panel 
of figure ||, whereas a higher-energy optical flux is greater 
than that of a lower-energy flux in the intrinsic spectra in 
the absence of lensing (see figure ||). This indicates that 
the amplification of a lower-energy flux has already been 
appreciable at t = —120 day. This results in the color 
change; i.e. chromatic feature is produced. 

3.4- Parameter dependence 

Figure [?] shows the light curves for the case with a 
smaller lensing mass by one order and the same physical 
impact parameter (We) (Case 2). Namely, (b,M\ ens ) = 
[0.01 (M /M lens ) 1/2 , O.IMq]. We see that if the lens mass 
is small, the peak magnification is reduced but the vari- 
ability timescale is near the same at this range, while the 



Fig. 7. Same as figure |6| but for the case with smaller lens- 
ing mass; (b, M lons ) = [O.O^Mq/M^) 1 / 2 , O.IMq)] (Case 
2). 

chromaticity of the standard disk still remains. This is 
more clearly shown in table 1. 

Next, in figure ^|, we assign a larger impact parameter; 
(b, Mi ons ) = (0.1, I.OMq) (Case 3). In this case, not only 
the total amplification become small but also chromatic- 
ity of standard disks no longer holds. Consequently, when 
we try to discriminate the disk structure, our success will 
depend not on the lens mass, but on whether the impact 
parameter is small or not. We will discuss this issue later 
(in section 4). 

In table we summarize characteristic values of these 
three light curves. Achromatic and chromatic features are 
clear in this table. 

Finally, figure ^ shows the dependence of the spectrum 
on the lens mass. The smaller the lens mass is, the nar- 
rower the resulting shapes of the spectrum for a given flux 
value become. This is because when the lens mass is small, 

1 /2 

the Einstein-ring radius is small (9e oc M lc ' ns , see Eq. [5]), 
so is the size of the region undergoing microlensing ampli- 
fication. This reduces the total light amplification. 

4. Discussion 

We have demonstrated how spectral behavior and mi- 
crolensing light curves depend on the disk structure. We 
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Case 


disk model 


band 


l°g -f(200d) 


l°g -f(l()0d) 


log F (od ) 




log[-F(iood)/-F(20od)] 


l°g[-f(()d)/-F(l()0d)] 


Case 1 


standard 


U 


-26.179 


-25.883 


-25.174 




0.296 


0.709 






K 


-26.288 


-26.061 


-25.766 




0.227 


0.305 




ADAF 


U 


-33.113 


-32.820 


-32.139 




0.293 


0.681 






K 


-32.885 


-32.592 


-31.910 




0.293 


0.682 


Case 2 


standard 


U 


-26.667 


-26.380 


-25.674 




0.287 


0.706 






K 


-26.775 


-26.555 


-26.264 




0.220 


0.291 




ADAF 


U 


-33.600 


-33.317 


-32.638 




0.283 


0.679 






K 


-33.372 


-33.089 


-32.410 




0.283 


0.679 


Case 3 


standard 


U 


-26.340 


-26.247 


-26.203 




0.093 


0.044 






K 


-26.446 


-26.352 


-26.311 




0.094 


0.041 




ADAF 


U 


-33.273 


-33.180 


-33.137 




0.093 


0.043 






K 


-33.045 


-32.952 


-32.909 




0.093 


0.043 



Table 1. The calculated fluxes at 3650A ([/-band) and 22000A (JT-band) at different epochs for case 1,2, and 3. 
Second, and third columns indicate the adopted accretion disk models and the energy band. Fourth, fifth, and sixth 
columns show the observed fluxes at t — 200 d (or —200 d), 100 d (or —100 d), and d. Seventh and eighth columns 
display changes in the fluxes from t = —200 d to t = —100 d, and those from t = —100 d to t = d, respectively. 
These two values are good indicators of the chromatic feature. 
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Fig. 8. Same as figure || but for the case with larger impact 
parameter; (b,M lens ) = (0.1,1.0M Q ) (Case 3). 



considered two representative disk models: the standard- 
type disk emitting predominantly optical-UV photons, 
and the optically thin, advection-dominated flow (ADAF) 
producing a wide range of photons from radio to X-7 rays. 
Using the microlensing events, we can distinguish the dif- 
ferent emissivity distribution of AGN accretion disks and 



will be able to map them in details and reconstruct them 
from the light curves. 

At the peak of the microlensing event, when the inner 
part of the accretion disk is largely amplified, distinctions 
between the two disk models are most evident (especially 
when the impact parameter is small). Thus, broad band 
photometry will be able to detect the color changes as 
shown in figure |], thereby revealing the structure of AGN 
accretion disks. 

Fortunately, such observations do not require good 
time resolution; for example, once per night like the MA- 
CHO project (e.g., Alcock et al. 1998) is sufficient, as 
long as a lens is located far from the observer. A typi- 
cal timescale of microlensing event is characterized by the 
Einstein ring crossing time of a lens object expressed as 



vt/D. 



6 xl0 b 



ol 



M, 



lens 



Mr. 



1/2 



200km s" 



7 M 

y Wkpc J 



1/2 



s. 



(15) 



where D means D\ S D Q \/ D os . For MACHO-like events 
when all the factors in the parentheses of eq. (JlJ) are 
of the order of unity, the event timescales become ~ 0.2 
yr. For the cosmological quasar case, on the other hand, 
we find roughly D ~ lGpc = 10 5 x lOkpc, and thus, the 
third parenthesis became large by factor ~ 300, yielding 
the event timescales being up to ~ 60 yr. Thus, in prin- 
ciple, a few years observation is sufficient to obtain the 
information about central disk structure. Practically, re- 
construction from microlensing light curve that is so-called 
'inverse problem' is not so easy, because, this requires a 
lot more to be done, e.g., take some (noisy, not perfectly 
sampled) light curve in a few filters. The source profile de- 
convolved by the observed light curve may comprise large 
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Fig. 9. Dependence of the microlensed spectra of the stan- 
dard disk (upper) and the optically thin flow (lower) on 
the lensing mass, Mi cns = O.IMq (dotted line), l.OAf© 
(dashed line), and 10.0M Q (long dashed line). The impact 
parameters are O.O1(M0/Mi ens ) 1 / 2 . The spectrum without 
microlensing is also plotted by the solid line. 



errors, if the light curve sampling interval is much longer 
than the event duration timescale, or if the accuracy of the 
observed magnitude (i.e., flux from the source) is much 
poorer, as stated before by Grieger, Kayser, & Schramm 
(1991). Therefore, we need good time resolution not for 
discriminating accretion disk model, but for mapping de- 
tailed accretion disk. This reconstruction of accretion disk 
structure may be done in our future work. 

To summarize the distinct properties of the two types 
of disks, we plot in figure |h] differences in the amplification 
factors at two different frequencies. Open and filled circles 
represent the cases with b = 0.1 and b = 0.01, respec- 
tively. The lensing mass is assumed to be Miens = 1-OM . 
As the angular separation between the lens and the disk 
center (it) decreases, AF increases (thus the point moves 
in the upper right direction in figure 0). Obviously, the 
color changes are more pronounced in the right panel (be- 
tween the U- and K-bands) than in the left panel (between 
the U- and V-bands), and for smaller impact parameters. 
We can easily discriminate the disk structure by compar- 
ing the optical and IR fluxes if b ^ 0.03. Furthermore, as 



clearly seen in figure |10|, the wider the separation between 
two observed wavelengths (or frequencies) are, the more 
apparent differences between the light curves of the two 
disk models become. Owing to this fact, the wider bands 
are preferable for more effective disk mapping. 

As a first step to understand the microlensing phe- 
nomenon, we have considered microlensing by a single 
star in the present study. However, the microlensing events 
caused by a single lens, especially the case in which dif- 
ferences between the light curves of the two disk models 
are apparent (i.e., b < 0.1), seems not so frequent for the 
actual objects. We thus need to evaluate the average time 
interval (At) between microlensing events with an impact 
parameter smaller than b. From Paczyhski (1986), we find 
that (At) is proportional to If we regard the event 
that the impact parameter is b < 1.0 as the microlensing 
event, one microlensing event which we are interested in 
will occur among 10 (b = 0.1) to 100 (b = 0.01) events. 

For example, we, here, assume the case of 
Q 2237+0305, although, in this case, single-lens as- 
sumption is not very good for this object because of 
its relatively large optical depth for microlensing (e.g., 
T = 0.4 ~ 0.7 from Rix, Schneider, & Bahcall 1992). 
It is widely believed that microlensing does occur in 
Q 2237+0305 at a rate of at least one event per year, i.e. 
(At) ~ lyr (Corrigan et al. 1991; Houde & Racine 1994). 
Hence, the microlensing event which we would like to 
observe is expected to occur once per 10 yr (for b = 0.1) 
to 100 yr (for b = 0.01). 

To conclude, microlensing events with such a small im- 
pact parameter may not be frequent, but once the mi- 
crolensing with a small impact parameter occurs, a large 
light amplification is expected. We have demonstrated 
that an optically thick accretion disk exhibits more cu- 
rious features in its spectral shape than an optically thin 
accretion disk. Differences between these two disk models 
are, essentially, caused by different energy balance equa- 
tions. In an optically thick accretion disk, generally, vis- 
cosity is small, accreting gas stays in the accretion disk 
for a much longer time than the disk rotation periods, so 
that gravitational potential energy of accreting gas can 
be effectively transformed into radiation on its way to the 
central compact object (in the present case, a supermas- 
sive black hole) . In an optically thin accretion disk, on the 
other hands, since viscosity is high, accreting gas falls into 
the central compact object rapidly in timescale on the or- 
der of the rotation period, and gravitational potential en- 
ergy of accreting gas is not so effectively transformed into 
radiation (e.g., Kato, Fukue, & Mineshige 1998). Thus, 
most of gravitational potential energy of accreting gas is 
advected into the central compact object and is not radi- 
ated away. These major discriminant physical processes in 
the accretion disks can be distinguished by our microlens 
diagnostic technique. 

In this paper, we assume that accretion disks are face 
on to the observer. In real situation, however, accretion 
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Fig. 10. Differences in the amplification factors at two different frequencies for the case of the standard disk (by the 
dotted line) and the optically thin flow (by the solid line), respectively. The adopted two frequencies are from the left, 
the U-band and the V-band (left panels), the U-band and the I-band (middle panels), and the U-band and the K-band 
(right panels), respectively. Open and filled circles represent the cases with b = 0.1 and b = 0.01, respectively. The 
lensing mass is assumed to be M lens = l.OM . 



disks may have non-zero inclination angles (i.e., i ^ 0). It 
will be complicated to quantitatively treat the dependence 
of inclination angles because of the increase in the number 
of parameters, such as the angle between the path of the 
lens object and the apparent semi-major axis caused by 
a non-zero inclination. Qualitatively, non-zero inclination 
affects timescale of the microlensing event in such a way 
to reduce the apparent disk size by factor ~ cosi. This 
effect seems to be small except the case that the disks are 
nearly edge on (i.e., i ~ 90°). This inclination effect will 
be systematically examined in a future work. 

Finally, recent high-resolution observations, including 
those from HST, have revealed circumnuclear starburst re- 
gions in AGNs (see Umemura, Fukue & Mineshige 1997, 
and references therein). If starburst regions exist around 
an accretion disk, microlensing event of AGN will make 
an interesting light curve; large-amplitude, short timescale 
variation (on a few months) caused by the microlensing 
of the accretion disk will be superposed on more gradual 
light variation (over a year) caused by the microlensing of 
starburst region (see figure pj). As shown in figure || the 
size of circumnuclear starburst regions is a tenth or hun- 
dredth as large as the central accretion disk. Hence, the 
microlensing timescale of circumnuclear starburst region 
is one tenth or hundredth as long as that of the accre- 
tion disk and the total amplification of the circumnuclear 
starburst region is smaller than that of the accretion disk 
from a rough estimate (see eq (2.28) of Schneider, Ehlers, 
& Falco (1992)). Thus, contributions of these two compo- 
sitions to the microlensing light curve will be able to dis- 
criminate, and long-term monitoring may provide a possi- 
ble tool to elucidate the extension of circumnuclear star- 
bursts although the presence of intrinsic variability tends 
to make such analyses complex. 




Contribution of 
Accretion Disk 




Contribution of 
Starburst Region 



• months 
- years 



— > 
Time 



Fig. 11. Schematic view of a microlensing event of an 
AGN disk and circumnuclear starburst regions. Upper fig- 
ure shows the situation of microlensing of accretion disk 
plus its circumnuclear starburst region. Lower figure shows 
its expected light curve. 



5. Summary 

The single-lens model as adopted in this paper is a good 
approximation, unless the apparent separations among 
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lens stars are appreciably smaller than the radius of the 
Einstein ring in the source plane. But, as we stated above, 
we should note that the images of Q2237+0305 could be 
subject to microlcnsing by multiple stars of the foreground 
galaxy. That is, we must consider, as a next step, the mi- 
crolensing events caused by 'caustics' which is produced 
by multiple microlensing. Such events may occur more 
frequently with shorter durations and steeper magnifica- 
tion (Wambsganss & Kundic 1995). Light curves by the 
caustics will show even more interesting features than the 
present case (Schneider & Weiss 1986). More realistic, ana- 
lytically approximated 'caustic' calculation has been done 
for the case of Q 2237+0305 (Yonehara et al. 1998). 

This diagnostic method for the central region of 
AGN can be used in other sources (e.g., MG0414+0534, 
PG1115+080, and the so-called 'clover leaf, H1413+117 
etc.). Future systematic and wide survey of quasars will 
discover other multiple quasars which our microlens diag- 
nostic technique can apply. 
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